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Abstract 

> 

We show how the fully resummed thermal pressure is rendered ultraviolet finite by standard 
zero-temperature renormalisation. The analysis is developed in a 6-dimensional scalar model 
' that mimics QED and has N flavours. The N — > oo limit of the model can be calculated 

completely. At a critical temperature, one of the degrees of freedom has vanishing screening 
mass like the transverse gauge bosons in four-dimensional finite-temperature perturbation 
theory. The renormalised nonperturbative interaction pressure of this model is evaluated 
numerically. 



OO 

a: 

1 Introduction 

The perturbation series for the pressure in finite-temperature QCD suffers from severe infrared prob- 
lems. In principle, these may be cured by a resummation technique. This resummation is most simply 
carried outW before renormalisation. It goes without saying that, after renormalisation, the pressure 
must be finite if it is to make physical sense, but it is far from obvious how the mathematics takes care 
of this. In this paper we study this in a model that may be regarded as a simplified mimic of QED, 
in which there are N particles each having the same mass and "charge". Our analysis makes use of 
standard techniques for renormalising composite operators. We find that indeed the usual renormal- 
isation, carried out purely at the zero-temperature level and therefore introducing no new quantities 
needing to be determined by experiment, renders the pressure finite. 

For simplicity, our discussion begins with the N — > oo limit of the model, which can be calculated 
exactly. The model is richer than the large- iV </> 4 theory which we have studied previously^ , in that 
now the self-energies in the large- iV limit vary with momentum, and wave-function renormalisation is 
needed. After renormalising the expression for the pressure and going some way towards evaluating 
it analytically, we complete the calculation of the large- N case numerically. 

Such a calculation is potentially useful even when N is not large. The formula for the resummed 
pressure involves the thermal self-energies of the fields, which inevitably are calculated in some ap- 
proximation from a finite number of Feynman graphs. The resummation then effectively converts this 
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finite set of graphs to a contribution to the pressure from an infinite number of graphs. While the 
exact form of the pressure must be ultraviolet finite, it is not obvious that it is still finite when only 
a partial set of graphs is included in the self-energy. One way of selecting a consistent approximation 
is to use in the finite- N case only the set of graphs that would survive to some given order in iV _1 if 
one were to take the large- N limit. In a theory in which there are a number of real unrenormalised 
fields 4> r , the pressure at temperature T is calculated' 1 ] from the thermal averages of the composite 
operators 4> 2 r : 

0^rP{T) = -k((M*)M0)) T ~ (o\Mx)Mo)\o)) x=o 



fl.la) 



'Or 

or, equivalently, 

__P(r) = -l(<T^ T .(x)^(O)> T -(O|T^(x)0 r .(O)|O)) ;c=o (1.16) 

Here, the differentiation is with respect to the unrenormalised mass of the field <p r , with all the other 
unrenormalised parameters kept fixed. One way to integrate this to give P(T) is to write 

m'^. = xml r r = 1, 2, . . . 

^(T)=y>^P(X) (1.2) 

Then integrate with respect to x from 1 to oo and insert the boundary condition that the pressure 
should vanish when x = oo, that is when all the masses are infinite*. However if, as is the case in 
QED, it happens that taking just one of the masses — mo s say — to infinity switches off all the 
interaction, there is a simpler method: 

P{T) = - / dm g£^P{T)+P {T) (1.3) 
Jml s om 0s 

where Po(T) is the contribution to the thermal pressure from all the fields except (f> s , with the inter- 
action between them switched off. 

The thermal averages of the composite operators that appear in (1.1) may be expressed as integrals 
over thermal Green's functions: 

(Mx)Mo)) T \ x=o = J $rD%(q) 

(TM*)M0)) t I x=Q = J (£f^(9) (1.4) 

The notation D^j,, is that of real-time thermal field theory in the Keldysh formalism W with a 
time path with a = 0: they are elements of the matrix propagator 



D rT (o) = M T (q°) 

where 

D rT (q) 



D rT {q) 
D* rT (q) 



M T (q°) (1.5a) 



g2- m 2 r _n rT (gO,q2) 



Our methods resemble those of a renormalisation-flow analysis' 3 '. 
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Here, Tl r T is defined in terms of the thermal self-energy matrix 



y>/ T 



e 

e \\i°\/ T 



rT 



m rT o 
o (-iu rT ( q )y 



M-\q^ 



Inserting (1.5) into (1.1a) gives 
d 



dm\ r 



1 ^ 2 J (2vr)"l (? 2_ m 2 r _ nrr(g )q 2 



ml-Il r ( q i)} 



(1.56) 



(1.5c) 



(1.1c) 



where H r (q 2 ) is the zero-temperature self energy and n(q) is the Bose distribution (e' 9 °'/ T — 1) 1 . The 
version (1.1b) of our basic formula gives instead 



d 



dmf jr 



P m = i f { im ^ 

1 j 2 7 (27r)«l q 2 -ml-U rT (Q ^ 2 



+ 



nir 



v 0r ~ n r (? 2 J H — '"-Or 

The equivalence of (1.1c) and (l.ld) may be seen from the fact* that each of the last two terms in 
the integrand of (l.ld) is the analytic continuation in q° of the real-time propagator, and by making 
a Wick rotation, so verifying that the integral over each of the last two terms is real. 

We apply the formula (1.1) to a mock electron-photon interaction in which, for simplicity, the fields 
are scalar; its unrenormalised form is 



nir 



U rT (q 



°,q 2 )} 



(l.ld) 



-INT 



(1.6) 



The masses are tooi for the electrons, and mo2 for the photon. In proper QED, C-parity or spin 
conservation removes one-photon-reducible graphs from the pressure; here we achieve this instead by 
making the photon an isovector, and the electrons isodoublets, so that a runs over 3 values with 
T a the Pauli matrices. We take space-time to be 6-dimensional, so that Ao is dimensionless and the 
divergences are similar to those of proper QED. In intermediate steps we use dimensional regularisation 
with n = 6 — 2e. 



There are iV identical electrons: the index r runs over N values. In the next section we set 

9o 



A = 



(1.7) 



and consider the case of large TV. The more general case, where TV is not necessarily large, is the 
subject of section 3. In section 4 we return to the large-N version of the model and evaluate the 
leading term in the interaction pressure completely. The thermal "photon" spectrum turns out to 
involve negative corrections to the mass such that there is a critical temperature, where screening 
disappears while keeping the plasmon mass nonzero. Right at the critical temperature, where the 
nonperturbative interaction pressure is still well-defined, the spectrum of our model is even rather 
similar to that of perturbative four-dimensional gauge theories in that it has a vanishing screening 
mass like the magnetostatic modes. Section 5 is a summary and discussion. 



* Formula (2.68) of reference [3] 
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(a) (b) 
Figure 1: self-energy graphs 



2 Large N 

In the large- N limit, the free- field pressure is linear in N, and the correction from the interaction is 
of order N°. To calculate this we need the leading terms in the photon self-energy g 2 ,7r5 ab and the 
electron self-energy, which are respectively of order N° and iV -1 and correspond to the graphs of 
figure 1. In figure lb, the photon propagator is the Dyson-resummed propagator with the photon self 
energy g^ir of figure la. 

If the photon mass mo2 — ► oo the interaction is switched off, and so we can use the version (1.3) of 
the formula for the pressure. With (1.1a) and (1.5) this reads 



P(T) =P (1) (T,mg 1 ) - \C Im [°° drug [ 

Jm 02 J 



(27T)6-2e\ g 2_ m g2 2 _^ 7rT ( g )q 2) 



1 



(2.1) 



q 2 -mg -g%TT(q 2 )) 

where C = |tr r 2 = 3 is the number of "photon" field components in our model. The integration over 

the photon mass m' 02 is with the coupling g and the electron mass moi fixed. P (1) (T, moi) is the 
contribution to the pressure from the electrons when the interaction is switched off. The N doublet 
fields tp r correspond to AN real fields, thus: 



P«(Tm 2 )-- 27V Im r dm' 2 I * q f l±Mg) ^1 

=AN7T I w^ n{q)e{q2 



(2.2) 



2 2 n 



Here, and elsewhere if needed, we assume the usual ie prescription. 

The next step is to express (2.1) and (2.2) in terms of renormalised quantitities. We introduce a 
renormalised photon mass m' 2 for each value of the bare mass m' 02 . The most appropriate scheme is 
on-shell renormalisation, in which the renormalised photon mass is given by 

™'i = m'g + glTr{m' 2 2 ) (2.3) 

The electron mass m 01 is renormalised similarly, though its renormalisation vanishes as iV -1 when 
N -> oo. 

We are going to calculate the pressure up to terms of order iV , so that in the calculation of Tr(q 2 ) and 
TTxiq , Q 2 ) we nee d make no self-energy insertions in the internal electron lines; these functions depend 
just on q and on the unrenormalised electron mass moi, which we may equate to its renormalised 
value mi. However, the renormalisation of moi will be important below. To leading order in N, the 
renormalised coupling is 

f{m' 2 ) = Z 2 {m' 2 )g 2 

l-^Tr'(m^) dm' 2 2 
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where 7t'(to 2 ) denotes dit{q 2 )/dq 2 evaluated at q 2 = m 2 . Thus, for fixed g%, 



-2 a 9 2 

9 ( ™ 2 ^ = 1 + g 2 (ir'(m 2 ) -ir'{m' 2 )) 

g 2 =f(m 2 ) (2.46) 

The equations (2.4) give pathologies reminiscent ^ of those of large- N (f> 4 theory. Similar pathologies 
were discussed for an exactly solvable model a long time ago' 5 ] . If we insist that g 2 . > 0, then we find 
from (2.4a) that g 2 — > as e — ► 0. However, if we take the view that the value of g 2 is irrelevant 
for physics, and simply choose g 2 to have some positive value, we find that g 2 (m 2 2 ) has a Landau 
pole. But provided that g 2 <C 1927T 3 , the pole is at such a large value M 2 of m' 2 2 that it is physically 
irrelevant: M 2 ~ to 2 exp(1927r 3 /g 2 ). Then our thermal field theory is a sensible theory provided we 
restrict ourselves to temperatures T <C M. 

The curly bracket in (2.1) is 

1 + 2n(q) 1 I 



■>«){■ 



(2.5) 



. q 2 — m' 2 — g 2 (m' 2 2 )iTT((f' , q 2 , m' 2 ) q 2 — m' 2 — g 2 (m 2 )Tr{q 2 ,m 2 ) / 
Here, ttt and w are convergent functions in 6 dimensions, when e — > 0: 

vfT(g°,q 2 ,m 2 2 ) =vr r ( g °,q 2 ) - {q 2 - m' 2 )*' (m 12 ) - ir(m> 2 ) 

7r(q 2 ,m' 2 ) = 7r(q 2 ) - (q 2 - m' 2 )7r'(m' 2 ) - ir(m' 2 ) (2.6) 

Of course, both ttt and 7f depend also on to 2 , but this is kept fixed at its physical value in both these 
functions. We may use (2.4) to change the integration variable from TOq 2 to m' 2 . Because we are 
taking g 2 < and g 2 > 0, when TOq 2 increases from TOq 2 to oo we find that m' 2 decreases from m 2 to 
— oo. These negative squared masses are just a calculational device and they do not enter in the final 
result but, provided that m 2 < 2toi, they mean that the squared renormalised mass and coupling 
defined in (2.3) and (2.4) are real throughout the integration, and the integral in (2.1) is 



_i rT r°° , a f d?~ 2£ q f i + Mg) 

2 Ll 2 J (27r)6- 2 4 9 2 -TO 2 2 -3 2 (m 2 2 )7f T (^,q 2 ,TO 2 2 ) 

1 



(2.7) 



q 2 — to' 2 2 — g 2 (m' 2 )ir (q 2 , to' 2 2 ) 

We may perform the mass integration, because from (2.4) 





dm' 2 



g 2 { m ' 2 )=g\m ,2 )K"{m' 2 ) (2.8a) 



and from (2.6) 



' 7f T (a°,q 2 ,TO 2 2 ) = -^vf( 2 ,TO 2 2 ) = -(g 2 - m 2 2 )7r"(TO 2 2 ) (2.86) 



dm'2 ' ' dm' 2 2 



From (2.8a), (2.8b) we find that (2.7) is 

d 6-2e, 



/d q 
j^z^n{q)6{q 2 - m 2 ' 
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d 



■6-2e, 



-^ CIm ./ (2^r (9)1 ° g m^q~ 2 +1 ° g g H(q 2 ,m 2 )+m 2 -q 2 \ (2 - 9) 

This integral diverges when e — > 0. The divergence must be cancelled by a similar one in Pq 1 ^ (T, mgi), 
given in (2.2), which we rewrite as 

P (1) (T,m2 1 )=p( 1 )(T,m 2 ) + J B(T) 

/j6-2e£, 
j^y-^n(k) [6(k 2 - ml) - 9(k 2 - O] (2.10) 

Evidently Pq^T, mf) is convergent, but B(T) is not. Because m 2 — goes to zero as iV _1 as A?" 
becomes large, 

/d 6 ~ 2e k 
^—^n(k)2^{k 2 -ml) (2.11a) 

Figure lb gives 

(m\ - m 2 01 )N = « 5 2 / ( £ )6 l g2 _^_^ (g2>m2) (fc _ g )2_ m 2 ! fe2=m? ( 2 - 116 ) 

So 

/^6-2e i 
(2 J )°-Vml- 9 wX) ' r(,W) <2 ' llc) 

/ f ]6-2el, i 

(We have used the fact that B(T) is real). 

We must show that the divergent parts of (2.9) and (2.11c) cancel when e — ► 0. From the graph of 
figure la, 

<f) = 2,/ (1)6-1 k 2 l m 2 {q _ k y _ m? ( 2 - 12 «) 

Also 

4 1 (g°,q 2 ,^l)-^(g 2 ,m|)= 2 6 _ 2e | d 6 " 2 ^ d 6 " 2 ^ ^(fci+fe-g) 

f77.(fci)27r(5(A; 2 - m 2 ) n(k 2 )2Tr5(kl - ml) Xo r/l2 2\n r/i2 2\1 /« 10 n 

1 , 9 9 — + , 9 9 — - tn(k 1 )n(k 2 )2Tr5(k( - ml)2Tr5(k 2 - m{) \ (2.126) 

I k 2 - m{ k{-m\ > 

and from this we may calculate ttt(q°, q 2 , m 2 ) because (1.5c) tells us that 

Re n rT = Re 

From (2. lid), (2.12b) and (2.13) it is immediate that 

Re 7r T (<7°,q 2 ) = Re {vfr(g°, q 2 , m|) - vf(g 2 ,m 2 )} (2.14a) 
It is familiar t 6 ' that Im ir(q 2 ) is even in q° and 

Im 7r(g 2 ) = -9{q 2 - Ami) ^ 6 -2e f d 6 - 2e k 1 d 6 - 2e k 2 5(k 1 + k 2 -q)2TT5(k 2 -ml)27r5(k 2 -m 2 1 ) (2.15a) 

From (2.12b) we see that 



Im {7t 1 T \q°,q 2 ,m 2 2 )-7r(q 2 ,m 2 2 )} = 

" (2^27 / d 6 ' 2 '^ 6 ' 26 ^ 6 ^ + fc 2 - q)27T5(k 2 - m 2 )2n5{k 2 - m?){n(fci) + n(k 2 ) + 2n(fc!)n(A; 2 )} 

(2.156) 

This is nonzero for both q 2 > 4mf and q 2 < 0. For q 2 > Am\ the (5-functions demand that both 
and k 2 are positive or both negative, so that 

ra(fci) + n(fc 2 ) + 2n(fci)n(fc2) = {1 + 2n(g)}{l + n(fci) + n(A; 2 )} - 1 (2.16a) 

and we find that the relation (2.14a) is true also for the imaginary parts of the functions involved. 
But for q 2 < the 5-functions require k® and k 2 to have opposite signs, and instead 

n(fci) + ra(fc 2 ) + 2n(fci)n(A&) = e(g )e(fc°){l + 2n(q)}{n(k 2 ) - n(h)} (2.166) 

In consequence, 

Im vr T (a°, q 2 ) = Im {vr T (g°, q 2 , to 2 ) - vf (a 2 , m 2 )} - 0(-g 2 )<M<A q 2 ) 

Mq°,cL 2 ) = o(q )Mq°,<i 2 ) + o(-q )M-q ,<i 2 ) 

%°)Mg°,q 2 ) = 4%°) ^ 6 _ 2e J d^k^-^k^k^^- q)2^5{k 2 -m 2 )2^5{k 2 -m^ 

(2.146) 

So finally, when the free-field pressure is subtracted off, the interaction pressure is 



r d 6 ~ 2e o r 
P(rr = -iCIm J ^pi^lo 



ff 2 7r T (o°,q 2 ,m 2 ) + to 2 - q 2 
to 2 , — q 2 



g 2 TT T (q°,c l 2 ,m 2 ) +m 2 - q 2 2 ir T (q°, q 2 , m 2 ) - vr(g 2 , m 2 ) - i6(-q 2 )<f> T (q , q 2 )} 



g 2 7r(q 2 , to, 2 ,) + m 2 , — g 2 g 2 — to 2 , — g 2 n(q 2 , to 2 ,) 

(2.17) 

For e — ► 0, that is in 6 dimensions, the integrand in (2.17) is finite, so it remains to check that its 
high-g behaviour is such that the integral is finite. From (2.12a) 

^ 2 ) = - 7T ^r7F(-l + e) / dx [m 2 - q 2 x(l - x)] 1 " 6 (2.12c) 
( 4 tt) Jo 



so that 



7f(g 2 ,m 2 ) = - 2 T(-l + e) I dx{ [m 2 -g 2 x(l-x)] 1 " e - [m 2 -m^(l-x)] 1 " £ + 
( 47r J Jo 1 

(g 2 - m 2 )x(l - x)(l - e)[m 2 - m 2 x(l - x)]" e } (2.12d) 

When g 2 becomes large, 7t(q 2 ,m 2 ) contains terms of order g 2 and g 2_2e , while the integral shown in 
(2.12b) that gives the difference between it^ and n is only of order 1/g 2 , so that ttt (g° , q 2 , m 2 ) ~ 
7f (g 2 , to 2 ,) and 

lQg g 2 7fT(g°,q 2 ,m 2 ) +m 2 - g 2 _ ^ 2 7f T (g°, q 2 , to 2 ) - Tt(q 2 , to 2 ) , c*(l/g 8 ) (2.18) 
g 2 Tt(q 2 , to 2 ,) + ml — g 2 q 2 — m 2 — g 2 Tt(q 2 ,m 2 ) 



On the other hand one can see from (2.14b) that when |q| is large, whether or not q° also is large, 
0r(g°,q 2 ) is exponentially small. Hence when q 2 is large the last two terms in the integrand of (2.17) 
combine to make the integral over q convergent. 

The central result of this section is our formula (2.17) for the interaction pressure to order N°. 
The reader may worry about some of the steps in its derivation, for instance an integration over 
negative renormalised squared masses for the photon in (2.7). Thus in intermediate steps we considered 
tachyonic photons! In the next section we will discuss the general renormalisation problem for the 
pressure. This will lead us in the 1/N expansion to another derivation of (2.17) which avoids these 
problems. 

3 Renormalisability of the pressure 

We now return to the interaction (1.6) and show how it leads in 6 dimensions to a finite expression 
for the pressure even when we do not take the large-iV limit. The form (1.6) is designed to simulate 
QED without introducing the complications that arise from spin, and so our analysis follows closely 
the familiar renormalisation of zero-temperature QED, such as is described in the book of Bjorken 
and Drelll 7 ] . As we shall see, the task of expressing the derivatives of the pressure with respect to the 
masses in terms of renormalised propagator and vertex functions leads us to a problem of overlapping 
divergences. This turns out to be similar to the overlapping-divergence problem for the vacuum 
polarisation in QED, where, following Dyson, one introduces a certain electron-positron scattering 
kernel (see chapter 19 of [7]). We shall follow a similar road here. 

Our notation will be as follows. We use the labels a, (3, . . . and a, b, ... to denote components of isodou- 
blets and isotriplets, respectively, and r,s, . . . for flavour labels. We also write the unrenormalised fields 
of the "electrons" and "photons" together as 

<f>lr = i'r 

r 2 = A a (3.i) 

These are the unrenormalised fields. As in the last section, their unrenormalised masses are moi, 
i = 1,2. Where it does not cause confusion, we will not explicitly write the isospin and flavour labels. 
The renormalised fields will be labelled with an additional suffix r; their masses are mj, which are 
related to the unrenormalised masses by 

m 2 - ml - U t (q 2 ) | =0 i=l,2 (3.2a) 

Here, Hi(g 2 ) and ^(g 2 ) are the self energies with the various Kronecker deltas factored off. From 
them, we also construct the two wave-function renormalisation constants 



Zi 



dq 2 



(3.26) 



J q 2 =m' 2 i 



and hence the renormalised fields 

4>iR = Z; l,2 <t>i (3.2c) 

We are assuming that < m\ < 4m 2 , so that the renormalisations are real. The propagators before 
and after renormalisation are 

D t ( q 2 )=i[q 2 -m 2 -U l (q 2 )]- 1 .l 
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D lR (q 2 ) = Z- l DM 2 ) = *[q 2 ~ m 2 - ^(a 2 )]" 1 • 1 (3.2d) 

where the two functions 

fi^g 2 ) = ^[n,(g 2 ) - n,(m 2 )] - (Z, - 1)(q; 2 - m 2 ) i = 1, 2 (3.2e) 

are finite and have value zero and zero derivative at q 2 = m 2 . In (3. 2d) 1 is to be read as 5 a p8 rs for 
i = 1 and as 5 a b for i = 2. 

We introduce the one-particle- irreducible vertex function j(p' ,p) which couples a pair of electrons of 
momenta p, p' to a photon. Renormalise it and the coupling, so that 

A7h(p',p) = A ZiZ 2 1/2 7 (p',p). (3.3a) 
We choose to fix 7r(p',p), and so define the renormalised coupling A, by imposing the condition 

(3.36) 



7R 

where 



SP 



7R 



= 7r(p,p') evaluated for p =p = (p — p') = — M (3.3c) 



for some fixed mass M. This renormalisation absorbs, for example, a divergence from the one-loop 
triangle graph (which we could neglect in the large- N limit). 



Note that the coupling (1.6) is invariant under the following C-transformation 

C: <j> lr — e(</>\ r ) T (r = l,...,N), 



irJ 

b 2 — > ~4>2 



where 



1 
-1 



(3.3c?) 



(3.3e) 



This forbids a nonzero vertex function for three "photons" in our model, similarly to Furry's theorem 
in QED. By the usual power counting arguments (compare [6]), we see then that the mass, wave 
function and coupling renormalisations (3.2), (3.3a)-(3.3c) make the (perturbative) theory finite. 

We need to introduce a 2 — ► 2-body connected scattering amplitude T. This is a 2 x 2 matrix connecting 
the channels 

1: 0i + 4>! 

2: 4>2 + 4>2 (3.4) 
Each element of T is also a matrix in isospin and flavour space. We define T to be amputated - there 
are no single-particle poles in its external legs - and we exclude from it all terms which are one-particle 
reducible in the s-channel, though T\\ does have a t-channel photon pole and T\i and T21 have t and 
u channel electron poles. We also introduce the two-particle-irreducible kernel K associated with T. 
It has no s-channel two-particle intermediate states and is related to T by 

T = K + KPT = K + TPK (3.5a) 

Here the 2x2 matrix P is diagonal; one diagonal element Pn is the tensor product of two electron 
propagators, and the other P22 is half the tensor product of two photon propagators. The 12 element 
of the first matrix equation in (3.5a) is drawn in figure 2, together with the definition of the matrix 
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Figure 2: one of the equations (3.5a) and the definition of the propagator matrix P 



r = 




Figure 3: Definition of the matrix of Green's functions that couple the composite fields to the 
electron-positron and photon-photon channels (3.4) 

P. The factor ^ in the definition of P takes account of the symmetry of channel 2 under interchange 
of the two photons. Because of the relation (3. 2d) between the renormalised and unrenormalised 
propagators, 

p r = z 1 pz- 1 



Z 1 

z 2 



(3.56) 



On the other hand, an amputated scattering amplitude on renormalisation acquires a factor Z 1 ^ 2 for 
each external leg, so the renormalised versions of T and K are 



T R = ZTZ 



K R = ZKZ 



(3.5c) 



Hence after renormalisation 



Tr — Kr + KrPrTr — Kr + TrPrKr 



(3.5d) 



Note that T and K have skeleton expansions that express them in a unique way in terms of integrals 
over the unrenormalised propagators and vertex functions in (3.2) and (3.3), while Tr and Kr have 
similar expansions in terms of the renormalised functions. 

We now define the composite operators constructed from the electron and photon fields 

X i(x) = <A(x)(t>i{x) - (0|4(*)<Mx)|0> 

X2 (x) = i[^(x)-(0|^(x)|0)] (3.6) 

Here, a sum over flavour and isospin indices is implied. Define also the one-particle irreducible Green's 
functions that couple these to the channels 1 and 2 defined in (3.4); they form a 2 x 2 matrix T, shown 
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graphically in figure 3. Each element of the matrix is a function of the momenta p and p' of the 
final-state particles. In fact for our purposes we need only consider the case p = p' , so that T is a 
function just of the single variable p 2 . To zeroth order in perturbation theory, r = Tq which, apart 
from Kronecker deltas for flavour and isospin indices is just the unit matrix: 











>ab 



(3.7a) 



The complete T has flavour and isospin structure 



f ii(p 2 )S rs 5 al3 t 12 (p 2 )5 ab 
r 21 ( P 2 )5 rs 5 a/3 t 22 (p 2 )5, 



ab 



(3.76) 



and it satisfies a Dyson equation 



r = r + t pt 
= r + tpk 



(3.7c) 



where we have used (3.5a). We define renormalised composite fields by 



Xr = Z x l X 



(3.8a) 



where the matrix Z v is to be chosen; then 



r R = z^rz = zr^oZ + r R p R K 



R 



(3.86) 



The isospin and flavour structure of T R is the same as that of F, that is (3.7b) but with scalar functions 
Tijn(p 2 ). Using the usual power-counting argument and the skeleton expansion for K R we see that 
indeed T R can be made finite by suitable choice of Z x , in such a way that the divergence from the 
last term in (3.8b) is absorbed in each order of perturbation theory (compare for example chapter 
19 of [7]). We also see from (3.8b) that any Z x that makes the elements of T R finite at one (but 
not necessarily the same) point in momentum space is an acceptable choice for this renormalisation 
matrix. 

In order to choose Z x appropriately, we return to the relation (3.2a) between the renormalised and 
unrenormalised masses. Remembering that IT(g 2 ) depend also on and m^ 2 , we have 



dm 2 
dm 2 0i 



d 

dm 2 0i 3 



(I 2 



_d_ 



U 3 (q 2 



dm 2 
#=m* dm 2 0i 



(3.9a) 



This gives, with the definition (3.2b) of Z i: 



dm 2 
dm 2 0i j 



s | d 

%3 dml 



n,(g 2 



Oi 



(3.96) 



The differentiation with respect to is with fixed bare coupling Ao- When applied to any Feynman 
graph for Ilj, it gives a sum of terms in which each internal line of type i in turn is doubled. We 
recognise therefore that the quantity on the right-hand side of (3.9b) is just 







dm 2 H 



U j (q 2 )=T ij (q 



(3.9c) 



11 



evaluated at q 2 = m 2 . Hence, if we define 



dm 2 



(3.9d) 



we see from (3.5b) and (3.8b) that 



'13 



(3.10a) 



Therefore the elements of r R (g 2 ) are finite when the final-state particles q are on shell, that is at one 
point in momentum space. From the arguments given above it follows that they are also finite when 
the particles are off shell. 

We are going to calculate the derivative of the pressure with respect to the renormalised masses at 
fixed bare coupling. It is convenient to introduce a renormalisation prescription in which, unlike in 
(3.2) and (3.3), fixed bare coupling corresponds to fixed renormalised coupling while the renormalised 
masses vary. In this modified prescription we replace the bare masses moi in (3.2a) with fixed bare 
masses /j,Qi, though keep the same bare coupling Ao- Then (3.2b), (3.2c) and (3.3) lead to a modified 
renormalised coupling A. We claim that A is a finite function of A. To show this, differentiate the 
definition (3.3) of A with respect to the bare masses at fixed Aq: 



dX 
dm\- 



X0Z1Z2 



d-y(p',p)\ 



SP 



+ A 



d 



dm\- 



logZi + -A 



1. d 



2 dm 2 3 



log Z 2 



(3.11a) 



Now 



Ac^- = 7>,P) 



(3.116) 



is the unrenormalised amputated Green's function for the fields Xji 0i> ^2 at zero four momentum 
for the Xj-field. Its superficial degree of divergence is —2, so it is rendered finite by the renormalisation 



1 . ^ 



fe = l,2 



(3.11c) 



Differentiating the definition (3.2b) of the Zi gives 



d 



dm 2 - 



\o g z l = { 

fc=l,2 ^ 



d 2 Ui(q 2 



(dq 2 ? 



d_ 

dq 2 



dm 2 



(3.11d) 



Then (3.11a), together with (3.9d), gives 



am^q 2 



1ir(p',P)\sp + A { ; + T^rn R (g 2 

dX , , . d 



dm 2 '^'^ ' " i (8q 2 ) 2 

72r(p',p)|sp + A^-f 2i R (g 2 )| g 



_d_ 

Oq 2 ' 



a 



dm 2 



+ iA 



+ 5 A ^2 r i2R(9 2 
lx (d 2 U 2 (q 2 ) d 



q 2 =m\ 



I {dq 2 ? 



+ ^r 22R (^ 



(3.11e) 



This expresses, in terms of renormalised quantities only, how A changes when the renormalised masses 
are changed, keeping Ao fixed. Thus, as we have claimed above, two renormalised couplings for different 
renormalised masses but the same Ao differ by a finite amount. The definition (3.3a,b) of A expresses it 
as a function of Aq and the bare masses men, which in turn can be considered from (3.2a) as functions 
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Figure 4: The diagrammatic expansion (3.12a) for dP(T)/dm 2 )1 



of Ao and the renormalised masses m r Then from A considered as function of Ao and the masses nii 
and the corresponding definition of A as function of Ao and the masses [ii we get two equations from 
which, in principle, we may eliminate Ao and so express A as a function of A (though in practice this 
will be a nontrivial task). However, from its definition it is clear that A remains fixed when the masses 
nii vary in differentiations with fixed Ao- 

We now return to formula (1.1). We choose the version (1.1b) and write it in the form 
d x f d n q 



^tr [roi^-'A}^)] t = l,2 (3.12a) 



Here T is defined in (3.7a) and 

Af T {q) = Df T {q) - Dj(q 2 ) = A$(q) • 1 (3.126) 

is the difference of the unrenormalised thermal 11-propagator and the zero temperature propagator 
for the field 4>j U = 1> 2). The trace is with respect to the isospin and flavour indices. 

The renormalisation of (3.12a) leads to a problem of overlapping divergences. To see this we consider 
the diagrammatic expansion of dP(T) jdm\^ , for example, shown in figure 4. The second term of figure 
4 can be interpreted as a correction either to the propagator on the right or to the vertex function 
of the composite operator \i on the left. This is analogous to the case of the vacuum polarisation 
function in QED (compare chapter 19 of [7]), as will be our methods for a proper renormalisation of 
(3.12a). 

We express Aj^ in terms of A^ R , which is defined in a similar way, but with the renormalised fields 
replacing the bare fields so that 

A&r(<?) = Z^A}},. (3.126) 
We use the fact that T is a constant matrix, but can be expressed from (3.7c) as 

r = r - tpk (3.126a) 

where the individual terms on the right-hand side are momentum-dependent. Now we choose as 
momentum argument on the right-hand side just the integration variable q of (3.12a). With this and 
the definition (3.9d) of the matrix [Z x ]ij, we have 

f(T) = -£ / (^^{[^(r-rPKjZ)]^ 1 -^^} (3.12c) 



dm 



or 

d x ^ f d n q 



dm 2 ; 



P(T) = " E / (0 tr { [( r * " r R P R K R )] y 2 1 -' Ag*} (3.12d) 
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Figure 5: The expression on the right-hand side of (3.12d) 

In using the formula, we must remember that the partial derivative is with the bare coupling Ao fixed, 
so the appropriate renormalised coupling is A. 

The right-hand side of the expression (3.12d) is shown graphically in figure 5. The renormalised vertex, 
kernel and propagator functions Tr, Kr, Dr, A^ r are, of course, finite for n = 6. But we still have 
three loop integrals to do, which we call Li,L 2 and L 3 , as indicated in figure 5. 

Let us first determine the degree of divergence for n = 6 of the loop integrations Li and L 3 . To this 
end, we write an operator-product expansion: 

<Ms)4r(°) = Ci(x 2 )l + doi* 2 ) XjR(O) + • • • (3.13a) 

3 

As the XjR have the same dimension as the left-hand side of (3.13a), the small-distance behaviour of 
Cij(x 2 ) is (x 2 ) , up to possible logarithms. The thermal propagator Dj^ K (q) is the Fourier transform 
of the vacuum expectation value of (3.13a). Because Aj^ R (g) is the difference between the temperature 
T and temperature propagators, the first term on the right-hand side of (3.13a) does not contribute 
to it and the leading power behaviour for q — > oo is obtained from the second term on the right-hand 
side of (3.13a). Thus, at worst, 

A lTR (g) ~ (g 2 )" 3 for q oo. (3.136) 

Power counting shows that the L\ integration is logarithmically divergent. For the second term in 
figure 5, we choose to do the loop integration L 3 first. Superficially, it is convergent. It really does 
converge if there are no divergent subintegrations. Of course, no divergent subintegrations occur in 
and Atr- Thus we could get a divergent subintegration only' 8 ' if the loop L3 is closed in Kr 
directly on one single skeleton vertex. But this cannot happen because Kr contains no s-channel 
1-particle-reducible diagrams. Thus we conclude that the L 3 loop is convergent; it gives a high-g 2 
behaviour 1/g 2 to that loop. From this, we see that the L 2 integration diverges logarithmically. 

So both terms in figure 5 are logarithmically divergent. It remains to show that these divergences 
cancel, leaving a finite result. From the discussion above it is clear that we will have to consider only 
the part of A^R(g) proportional to (g 2 ) -3 for q — ► 00. The higher terms in the operator product 
expansion (3.13a) lead to convergent contributions in all loops Li, L 2 , L 3 . Thus, for the discussion 
of the convergence we can replace A^ R in (3.12d) and figure 5 by any expression having the same 
(g 2 ) -3 behaviour for g — ► 00. We choose the following 4-point Green's functions: 

J^k(q, P )= I d n x I d n z 1 I d^2e^e- ipZl e ipZ2 {(0|T(^ R (x)(Al R (0)4R(^i)^R(^))|0) 

J J J 1 (3.16) 

where p is an arbitrary fixed momentum. For simplicity, we continue not to write explicitly the 
isospin and flavour indices of the fields. Those attached to (j)^ and 4>\r are carried by Jak, while 
those attached to i^r and </>£ R are equal to each other and summed. This summation excludes s- 
channel one particle reducible diagrams from Juk- Inserting here the operator product expansion 
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1 p 

Figure 6: Diagrammatic representation (3.17) of Juk{l,p) 



(3.13a), we see that the leading term for q — ► oo is again given by the CV,(x 2 ). In the loop integrals 
£1,2,3 with the thermal propagators AjT R possible divergences can only be proportional to the two 
thermal expectation values of Xj-r(O) in the operator product expansion (3.13a), where j = 1, 2. Thus, 
in order to prove convergence, we choose two "trial" terms J4 corresponding to the index k = 1,2 
in (3.16) which have linearly independent contributions from the Xj'r(O) C? = 1)2) in the operator 
product expansion. The structure of J^k is as follows (compare (3.5c), (3.5d) and figure 6): 

J4ik(q,p) = trJ?- 1 [P R (p)] ik (27T) n [5(q - p) + 5 k2 5(q + p)} + 

1 . (3.17) 

r- l [P R (q)T R (q,p)P R (p)] lk 2 k - 1 } 

where tr p stands for the summation over the isospin and flavour quantum numbers connected with 
the p-lines in figure 6. 

Let us then insert J 4 in (3.12d): 

E / ?Sr tr 4 [rR - rRPRKR] ^ 2l " JJ ^} 

i 

= E7^rN{[rR-r R p R K R ]^ g) 

i V ; (3.18) 
tr p [p R (p) ifc (2vr)™ (S(q - p) + 5 k , 2 5(q + p)) + [P R (q)T n (q,p)P R (p)] jk 2 k ~ l ] } 

= tr p {[r R ( P )p R (p) - (r R p R K R p R )( P ) + (r R p R T R p R )( P ) 

- (r R P R K R P R T R P R )(p)] ifc 2 fc_1 } . 

Using now (3.5d) we see that the last three terms on the right-hand side (3.18) cancel, leaving us with 
the first term which is obviously finite, being the product of the renormalised quantities T R (p) and 

Pr(p)- 

This concludes our discussion of the renormalisation of the derivatives of the pressure with respect to 
the masses. We have shown that in (3.12d) all divergences cancel leaving us with an ultraviolet finite 
result. 

It remains to make the 1/N expansion of the general formula (3.12d), where only renormalised quan- 
tities appear. This is straightforward but turns out to be a lengthy calculation. We give some details 
in the Appendix. As can be seen from there, starting from (3.12d) the unrenormalised masses and 
coupling are never encountered any more and no integration over negative renormalised squared mass 
is needed. The differential of the pressure with respect to both squared masses is obtained from (A. 5), 
(A.6), (A.23) as: 

J2dm^P(T,mlml) = E dTO ? A { P o%X) + ^V,™!) + ^(T, m?, ml) 
j ° m 3 j dm 3 L (3.19) 

+ 0(1/N)} , 
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P (1) (T,m?)=4iVvr J ^-n(q)6(q 2 - m 2 ), 
pf (T,m 2 ) =C*J-£jLn( q )0tf-m>). 



(3.20) 



Pq^T, m 2 ) is the pressure from the free electrons of renormalised mass mi, Pq 2 ^(T, m 2 ) similarly for 
the photons and (0)pint( T) 

m 2 ,m 2 ) i s the interaction pressure to order N° as in eq. (2.17). (To be 
consistent with the notation of this section we write here (°)_P int for the expression (2.17).) The mass 
dependences are indicated explicitly, and we remind the reader that always the coupling constant g is 
kept fixed. 

To summarise: We have calculated the derivatives of the pressure dP(T, m\ , m 2 )/dm 2 to order N 
and N° for all masses of photons and electrons satisfying < m| < 4m 2 . Starting from one point 
(mf^m 2 .) in this physical region we can integrate along any path C\ in the m 2 — m 2 . plane (figure 7) 
running to infinite masses, but always staying in the physical region. 

Using then as boundary condition that at infinite masses the pressure vanishes, we get: 

P(T,m 2 ,m 2 ) = - I V dm' 2 -^P(T,m 2 ,m 2 ) 

J Ci j dm j (3.21) 

} (T,m 2 ) +P (2) (T,m 2 ) + ^P int (T,m 2 ,m 2 ) + 0(l/N) 

On the other hand, in the simpler calculation performed in section 2 we had to integrate along C2, 
that is at fixed m 2 = m 2 with m 2 2 running from m 2 to —00. Clearly, C 2 leaves the physical region. 
Nevertheless the result obtained in (2.17) is correct, as we have shown here how it can be obtained 
without going to unphysical values for the masses. 

Thus the final result for the pressure to order iV 1 and N° is exactly as described in section 2: To order 
iV 1 we get the pressure of free "electrons" and "positrons" of renormalised mass mi, to order iV we 
get the free pressure of "photons" with renormalised mass m-2 plus the interaction pressure (°)p(T) mt 
of (2.17). 

4 Evaluation of the large- N limit 

We now turn to the evaluation of the exact result (2.17) which we obtained in the limit N — > 00. Since 
(2.17) is UV finite, we can put e — > 0. All formulae in this section thus refer to n = 6 dimensions. 



_ p(i) 
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Equation (2.17) can now be evaluated by a number of nested numerical integrations which involve the 
following functions as building blocks. Firstly, the zero-temperature contributions to the "photon" 
self-energy: 



Re ir(q 2 ) = —^3^ J dx [m\ - q 2 x{\ - x)]\og 



m 2 — q 2 x{\ — x) 



m 2 — m 2 x(l — x) 



+ l(q 2 -m 2 2 , 



(4.1) 



Im ir(q 2 



1 f 1 
32^ J dX 



1 — x) — m 2 ] 9(q 2 x(l — x) — m\) 



1 



192vr 2 



^-4„;)fa 2 - 4 "-' 3/2 



(4.2) 



which can be expressed in terms of elementary functions, though we have done so only for the imaginary 
part. 

Secondly, the thermal contributions. The real part of the thermal self-energy is given by 



Re (MM-^)) = T ^J~dk^^^ 



+ 
4 



/ 2m k q - q 2 
V 2fc|q| 



/ 2io k q + q £ 
V 2fc|q| 



log 
log 



u k qo - fe|q| - \q 2 



uj k qo + k\q\ - \q 2 
uj k q + fc|q| + \q 2 



u k q ~ k\q\ + \q 2 



2q 2 



k\q\ 



(4.3) 



where uj k = \Jk 2 + m\. The imaginary part is most easily obtained from 7r^ 2 via 

I mvrT = _l ie (g )(l_ e ^o )7r 12 

and using f(x)f(y) = f(x + y)[l + f(x) + f(y)} for f(x) = l/(e x - 1), yielding 



dk 



k 3 n 



M [ 



+ 



u k \ 

- q 2 
2&|q| 

( 2uj k qo + q" 
V 2fc|q| 




In (2.17) we also need frr- Its real part is identical to (4.3), whereas its imaginary part differs from 
(4.4) by replacing all the sign functions e with 1. 

In the limit |q| —>■ 0, the integral in (4.4) can be done analytically with the result 

Im 7t T (q , 0) = Im 7r T (q°, 0) = (l + 2n{\q )) Im 7f(g 2 = q 2 ) (4.5) 



Another limiting case which can be solved is the high-temperature limit \qo\, |q|, m\ <C T, which gives 

rp2 



24vr 



1 



4 
q 2 



, » , »+ q 

2 q qo ~ q 



+ 0(T) 



(4.6) 
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Remarkably, this is basically the same function that appears in the longitudinal component of the 
polarization tensor of hot QED and QCDM, except that here it comes with a reversed over-all sign.* 



As a consequence, the spectrum as read from the analytic structure of the full thermal "photon" 
propagator l/(q 2 + m 2 . — g 2 TTT(q , q 2 )) is rather unusual. 

In the case of initially massless "photons", the full thermal propagator still has singularities at the 
light-cone, because (4.3) vanishes at q 2 = 0; our massless "photons" do not acquire thermal masses. 
However, there are nontrivial corrections to the residues of the poles at q 2 = according to 



Z(q 2 ) = lim (1-g 2 « T y 



1 + 



9 



47r 3 q< 



Jo 



dk 



k 2 n(uj k ) 



2k °^ u>k — k 



(4.7) 



< 1 



This vanishes as |q|/T — ► and also as m\/T — > 0. In the high-temperature limit it is very small for 
momenta that are not at least comparable with T in magnitude. So while no mass gap is generated, 
as the temperature increases thermal effects progressively remove the modes with larger and larger 
momenta. In the infinite-temperature limit the residue of the pole becomes zero, so that then there 
are no propagating plasmons at all. 

On the other hand, for non-zero or non-neglible "photon" mass, there are always simple poles in the 
photon propagator. In this case there are thermal mass corrections, but they are negative, towards 
lighter (but nonzero) effective masses. At the same time, the residues of the corresponding poles are 
diminished. 

At qo = 0, ttt normally gives the screening mass-squared for static fields, which in gauge theories 
is the Debye mass. While screening corresponds to poles in the propagator for imaginary values of 
the spatial momentum, in our model we have a pole at real spatial momentum if the temperature is 
larger than some critical temperature T crit . For T » ra 1 ,m 2 this pole is located at q 2 = g 2 T 2 /(24ir), 
according to (4.6). A similar behaviour has been found in the gravitational polarization tensor of 
ultrarelativistic plasmas when evaluated on a flat-space background There the value of |q| at the 
pole at qo = is identified with the so-called Jeans mass characterizing the gravitational instability 
of the plasma. In our case, such a pole seems to reflect the fact that the potential of our model is 
unbounded from below, so that when the "photon" mass is small enough, thermal fluctuations can 
lead to a run-away symmetry breaking without the need of tunneling. 

For non-zero qo, there are no poles at space-like momenta, because for those (4.6) has a large imaginary 
part proportional to T 2 corresponding to Landau damping. So there are no propagating thermal 
tachyons. 

To summarize, the spectrum of our model in the limit of large N as read from the thermal propagators 
is the following. The only thermal corrections occur in the photon spectrum, which for temperatures 
sufficiently small compared to the photon mass 771,2 consist of negative (momentum-dependent) correc- 
tions to 7772- The latter are largest at low momenta and tend to zero for very high momenta. The cor- 
responding dispersion law is depicted in figure 8a for g = 10 and m\ = m 2 = m = T < T crit ~ 1.236m. 
Up to the critical temperature, the static fields have finite screening length, which becomes infinite 
at r cr ;t, whereas the plasma frequency (the long- wavelength limit of the dynamical mass) remains 
nonzero. Right at the critical temperature, the spectrum is thus very similar to that of the transverse 
vector bosons in the high-temperature limit of 4-dimensional gauge theories: a vanishing screening 



* This abnormal sign in </>|-theories has previously been noted in reference^. 
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Figure 8: The location of poles in the "photon" propagator for g = 10 at T = to = mi = m 2 (a), 
which is below the critical temperature, and at T = 1.5m (b), which is above it. The two dashed 
lines mark the light cone and the zero-temperature mass hyperboloid q'q = q 2 + m\. The line 
between them is the location of the poles in the thermal propagator. In (b) the second full line 
below the light cone marks zeros of the real part of the inverse propagator, which correspond to 
poles of the propagator only at its intersections with q = because of large Landau damping 
for q > and q 2 < 0. 



mass together with a nonzero plasmon mass. For temperatures above T cr ; t , there is a pole at space- 
like momentum q = and q 2 = to 2 signalling a Jeans-type instability As shown in figure 8b for 
T = 1. 5m > T cr it, the real part of the inverse photon propagator has zeros for space- like momenta 
with q 2 < m 2 , which gives rise to a pole of the propagator where also the imaginary part vanishes, 
which is at q = 0. 

The presence of a pole at qo = and |q| = mj > causes the interaction pressure 



P(T) 



int 



he 



+ arg 



d 6 q 



2n(q) arg 



g 2 ir T (q ,q 2 ,ml) +m% - q 2 



(2tt) 6 L 

g 2 7f T (g°,q 2 ,m|) + to 2 , - q 2 
g 2 ir (q 2 , to 2 ,) + to 2 , — q 2 



m\ — q 2 



+ g 2 Im 



7r T (^°,q 2 ,TO|) 



q 2 to 2 , g 2 Tr(q 2 ,ml) 



(2.17) 



to become IR singular for all T > T crit . In the first term, for |g | smaller than some finite number, the 
real part of g 2 TTT + to 2 , — q 2 changes sign from positive to negative as |q| is changed from some large 
value to a sufficiently small one (see figure 8b). The imaginary part on the other hand is negative 
throughout except for \qo\ = 0, where it vanishes. The argument in the first term therefore approaches 
a step function — ir9(mj — |q|) as \qo\ — > 0. This causes the first term in (2.17) to diverge logarithmically 
in the IR when T > T crit . 

For finite photon mass and T < T crit , the pole at qo = and |q| = toj > is absent and (2.17') 
appears to be well-defined. 

Turning finally to the numerical evaluation of (2.17), we can distinguish three different regions of 
integration depending on the appearance of imaginary parts. 
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The inverse photon propagator g 2 itT + m 2 ~ Q 2 nas imaginary parts only for q 2 > Am 2 (pair creation) 
and for q 2 < (Landau damping). The quasi-particles described by the time- like poles of the photon 
propagator are therefore undamped and stable, provided m 2 < 2m\. We denote their position by 
wr(q 2 )- F° r ko| < ^r(q 2 ) the real part of g 2 7Pr + m\ — q 2 is positive and for \qo\ > u>r(q 2 ) it is 
negative. 

Correspondingly, we have: 

I — spacelike momenta, q 2 < 0: Tt is real in this region, but ttt has imaginary parts corresponding to 
Landau damping, so that 



P(T)T* = -C I j™ dq {(2n(q ) + l)^g[g 2 7t T (q°,CL 2 ,m 2 2 )+m 2 -q 2 



•|q| 

,2i 



2 Im 7r T (g°,q 2 ,m|) -i 
q 2 — ml ~~ 9 2 Tf(q 2 , m 2 ) J 



(4.8a) 



/I — timelike momenta below threshold, < q 2 < Am 2 : the first two terms in (2.17) contribute only 
for 1/q 2 + m 2 > |<7o| > o>r(q 2 ), whereas the last term has a pole in this range with unit residue thanks 
to our on-shell renormalization scheme. The integration over q can be carried out with the result 

r rf 5 n ( p\/q 2 +™2 _ i / 

P{T) »=- C (2^^ 21 ° g e^)-l +"A<f)-^ 2 + m 2 

J { } 1 (4.86) 

2 Re 7r T (g° = ^/q 2 + m 2 , q 2 , m|) N 

2 Vq 2 + m 2 



III — timelike momenta above threshold, q 2 > Am 2 : 

P(T)fh = -C [ ^ dg {2n( (7 o)(arg[ 2 7f r ( (7 ,q 2 ,m 2 )+m 2 - (7 2 ]+7r) 

+ arg[g 2 Tt T {q () ,^,rn 2 2 ) + m 2 2 -q 2 ] - arg [g 2 7t(q 2 , ml) + m\ - q 2 ] (4.8c) 
7r T (g°,q 2 ,m|) 



+ 5 Z Im 



q 2 — m 2 — g 2 Tr(q 2 , m?,) 



The numerical evaluation of these expressions is quite challenging because of large cancellations among 
the individual contributions so that rather high working precision is needed. Our results for the 
nonperturbative interaction pressure are given in figure 9 for m = m\ = rri2 and three different values 
of the coupling, = 1, v^lO, and 10, and these are compared with the strictly perturbative contribution 
proportional to g 2 , which also has to be computed numerically because we do not restrict ourselves 
to the high-temperature limit and so the full momentum dependence of the thermal self-energy enters 
there too. 

For temperatures T <C m there is rather little difference between the perturbative and the nonpertur- 
bative results. As the temperature increases, the latter grow bigger than the former until they abruptly 
end in a singularity at the critical temperature. For sufficiently small coupling and for mi ~ m 2 , the 
critical temperature above which (2.17) ceases to exist can be estimated from the high-temperature 
expansion of the thermal photon propagator following from (4.6): T crit ~ V2Airm 2 /g ~ 9m 2 /g. The 
actual values for g = 1, \/l0, 10 are T crit ~ 9.133, 3.164, 1.236 times m, respectively. At these temper- 
atures the thermal pressure ceases to exist, because there is sort of a phase transition to a run-away 
and therefore inexistent broken phase. 
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P(T)' mt /(Cg 2 T 6 ) 




Figure 9: Numerical evaluation of the interaction pressure (2.17) for three values of g in compar- 
ison with the leading perturbative contribution. 

At exactly T = T CT a, we have a situation which is closest to gauge theories in 4 dimensions, because 
there mj = 0, corresponding to a vanishing screening mass as is the case for the magnetostatic modes 
in perturbation theory. The interaction pressure is still well-defined and is given by the end-points 
of the various curves in figure 9. A conspicuous difference from the results of 4-dimensional gauge 
theories is that the interaction pressure is positive, which is related to the abnormal sign of all thermal 
mass corrections in our model. 

5 Conclusions 

We have pursued two purposes in our study of our six-dimensional scalar model. Firstly, we have 
investigated how precisely the nonperturbative formula for the thermal pressure proposed in reference 1 
is rendered finite by standard zero-temperature renormalisation. 

The renormalisation of the derivatives of the pressure with respect to the masses leads us to a problem 
of overlapping divergences which we solved in a manner analogous to Dyson's method for QED, 
introducing a certain 2-2 scattering kernel. We derived the renormalised Dyson-type equations, which 
turned out to be essential for our discussion of renormalisability. Finally, we had to invoke Wilson's 
operator product expansion. We think that we can draw the lesson from this that in more complicated 
theories like QCD things will not be simpler and one will again have to deal with overlapping divergence 
problems. On the other hand, having an expression for the pressure in terms of renormalised Green's 
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functions as given for our model in (3.12d) and knowing the essential equations these Green's functions 
must satisfy in order to have ultraviolet finiteness, may help to devise consistent approximation schemes 
leading to finite results in all orders. With these methods we should also be able to study explicitly 
the effects occuring when a particle — in our case the "photon" — becomes unstable. In the Green's 
functions for T=0 this amounts simply to a pole moving from the first to the second Riemann sheet 
and thus we expect that our result for the pressure, which is expressed in terms of these functions, 
should not change drastically. 

Besides these general aspects, we have investigated the large N limit of our model which diagrammat- 
ically is similar to QED in the limit of large flavour numbers. In this limit the leading contribution 
to the interaction pressure comes from a ring resummation of the photon polarization function, while 
the electron lines remain undressed. When N is not large, keeping only this contribution corresponds 
to what is known as random-phase approximation (RPA)! 11 ! in many-body physics. In contrast to 
the simpler ring resummation of the Debye screening mass! 12 ] , one has to deal with a resummation 
of a momentum-dependent quantity. In practice, however, one usually aims at an (improved) pertur- 
bative scheme and uses this resummation only as far as needed to extract the next-to-leading order 
term in the interaction pressure, which because of the infrared singularities in the usual series is 
nonanalytiJ 13 ' [ 14 1 in e 2 . Here we have found that the RPA can be interpreted as the leading term in 
a large- iV-expansion and we have retained the full nonperturbative information that it incorporates. 

In our six-dimensional scalar model we have in fact encountered rather drastic resummation effects, 
because this model has a critical temperature above which an instability similar to the gravitational 
Jeans instability occurs. So despite the diagrammatic similarity with QED, this theory is rather dif- 
ferent from it. But below the critical temperature we were able to obtain a nonperturbative expression 
for the interaction pressure, and evaluate it numerically. Right at the critical temperature, where the 
nonperturbative interaction pressure is still well-defined, the spectrum of our model is even rather 
similar to that of perturbative four-dimensional gauge theories in that it has a vanishing screening 
mass like the magnetostatic modes. 

The computation of the nonperturbative interaction pressure in more realistic theories such as ordinary 
QED in the limit of large flavour numbers would be technically not too different from what we have 
done here. Similar simplifications seem to be of interest even in QCD in the small- A~ c large-iV/ limit! 15 !. 
We plan to investigate those theories along the above lines in a separate work. 
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Appendix 

In this appendix we show how the 1/N expansion of the general formula (3.12d) leads us to the results 
of section 2. In the first step we perform all the traces implicit in (3.12d). For this we define: 

K jlR (q,p) = i2 2 -i- l (tT q tT p )(K jlR (q, P )), (Al) 
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Figure 10: Diagrammatic picture for the expression (A.l) 



where tr g (tr p ) means the trace over the internal indices of the lines where the momentum q (p) flows 
(figure 10). We also define the differential forms: 



(A2) 



which can be expressed through the self-energy functions using the renormalised version of (3.9c): 

Z 



UJ 



i P ) = { - drnf^ [p 2 - m 2 - n,(p 2 )] + [p 2 - m 2 - IL> 2 )] dm 2 8 



/ — lo g#} (A3) 



dmf Zj 



Here Zj are the wave function renormalisation constants for the theory with the coupling A, but masses 
fij (compare the paragraph following (3. lie)). We find then for the pressure from (3.12d) 



d n p 
(2vr)" 



[q 2 -m* -Ujtf)]- 2 K jin ( q ,p)A% R (p)} 

3,1 



(A4) 



This is the starting point for the 1/N expansion. For any (scalar) quantity F we write F = ^ r ( r )p 
with ( r ^F the term of order (1/N) r . We get then from (A. 4) for the term of order N 



Zydm 2 ^ ^P(T) = -2Nj «V( g ) . ^ A^q) 

= -2Ndm 2 J ^-n( q )27T5( q 2 -m 2 ). 



(A5) 



Therefore, to order N we get the pressure from the free electrons and positrons of mass mi, that is 
P (1) (T,m?) of (2.10): 

(" 1 )p(T) = P (1) (T,m 2 ). (A6) 

To order iV° we get from (A. 4) after some work, and using the fact that the leading contributions to 
K are single-particle exchanges: 



Ydm 2 ® (°)P(T) =-f-£*{2N 
^ 3 dm 2 - v ; J 2vr )™ I 







+ dm\ (1) n;(m 2 ) 



(0) Alk(?) 



+ 2Ndm\ (1) A£ R (g)-±C ^ Z 2 







n 2(^2) - m i) 



(0) A^)}, 



(A7) 



23 



where, with the notation of section 2, 



'ILtf) = ^cj^[( g - k) 2 - mi]- 1 [k 2 -ml- gH (k 2 , ml)] ~' 



(0) n 2 (, 2 ) = 5o M<? 2 ), 



(A8) 

(A9) 
(A10) 



The unrenormalised quantities IIi,2, 2^, <?o occurring in (A. 7) are to be understood as functions of 
the renormalised parameters. Of course the individual terms in (A. 7) contain divergences, which must 
cancel in the sum as we know from sections 2,3 and as we will again see explicitly below. 



Now we use the general sum rule 



dq 



A£r( 9 ) - A# R ( g ) 



12 



(All) 



which we already mentioned in connection with (l.ld), to express the integrals over A 11 in (A. 7) by 
integrals over A 12 which in turn is related to the imaginary part of the thermal propagators. In this 
way we obtain 







dm ■ 



(o) 



p(T) = m +m + V3, 



(A12) 



where the differential forms 1/1,2,3 are given by 



771 = 2N 



f d n c 
J (2*) 







E*»3^? (1)n ^ a ) 



n(q)2ir5(q 2 - m\) 



k 2 =m\ 



772 = 2Ndm\ 



d n q 



Im 



[ + 2n{q) m 1T (Q) 



-c 



(2tt)" 1 (q 2 - m 2 ) 2 
: d7T ( m 2) ,2] f d n q 



dm 2 , — g 2 — - — 7^— dm 2 

omf 



Im 



(q 2 — ml) 2 
1 + 2n(q) 



(i) 



(A13) 
(A14) 



(2ir) n lq 2 — ml — g 2 7TT{q,ml) 



q 2 - ml - g 2 n{q 2 ,ml) 



with 



WlW) = Wn^g 2 ) - Wll^m 2 ) - (q 2 - m 2 )Wu[(m 2 ) 



(A15) 
(A16) 



g 2 C f d n k 



+ 



N 



f d n k f 

J (2^r 1 . 



(q — k) 2 — m\ 



+ n(q - k)2ir5{{q - k) 2 - m\) 



[k 2 - ml- g 2 TT T (k,ml)] 



2M-1 



{q-k) 2 -m 2 [k'-ml-g^ik^ml)]- 1 



— 1 



in(k) 



(q — k) 2 — m\ 



+ 2tt5 {(q - k) 2 - m 2 ) \n{k) (9(q°)6{k° - q°) + 9(-q°)6(q - k )) 



+ n(q-k)(6{q )6{-k ) + 6(-q )6{k°))^2 1m [k 2 - m 2 2 - g 2 7t T {k,m 2 2 )] '} (A17) 
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Here Wlli((/ 2 ) is the renormalised electron self-energy function to order 1/N and ^Uiriq) is the 
corresponding thermal function. 

The forms rji and 7/3 have already a simple structure. The form 772 which arises from the 1/N term 
W Aj^ R of the electron propagator in (A. 7) is more difficult to handle. The strategy is to insert (A. 16) 
and (A. 17) in (A. 14), which leads to integrals over q and k, and then to perform first the ^-integration. 
In this way we get after some nontrivial calculations: 







(1) n 1 K)^ ? + (1) n' 1 K; 



f] 2 = 2N dm\ 

+ ^ 2 dm\ J^lui{(l + 2n(k))[k 2 - 

- [k 2 -ml- g^k^ml)]- 1 ^(k 2 )}. 
Putting now everything together, we arrive at: 



d n k 
(2vr)" 



n(k)2ir5(k 2 - m\) 



m\ — <7 2 7r<r(fc, m\ 







dmf 



7T T (k) 



(A.18) 







2 

j dm 2 



dm 2 -q 2 CIm 
1 2^ 



tr) n \ l 



2 2 2-/2 2\T — 1 ^ 

q -m 2 -g 7r(q ,m 2 )J ^-j 
d 



(7r(q 2 ) - 7r(m 2 )) 



- [1 + 2n(q)] [q 2 - m 2 - g 2 Tt T (q, m 2 )] ^-^ (7r T (q) - Tr(m 2 )) 
dm^Clm J -^^{[q 2 ~m 2 2 - g 2 Tt(q 2 ,m 2 2 )]~ 1 
[l + 2n(q)\ [q 2 - m 2 - g 2 n T {q,m 2 2 )]~ 1 } 



In the first term on the right hand side of (A. 18a) we can use: 



(A18a) 



^itK) 

2Cg 2 Im ■ 
= 2Cg 2 Im 



dm\ 2 
d n q 1 



dm\ 



+ c.c. 



(27r) n (k — q) 2 — m\q 2 — m\- g 2 ir(q 2 ,m 2 ) 
d n q 1 f d n k 



k 2 =m 2 J (ZTT) 



(27T^q 2 -m 2 -g 2 7t(q 2 ,m 2 )J (2vr) 



n{k)27rS(k 2 - m() 



{k — q) 2 — m\ 



-Cg 2 Im 



d n q 



TT T {q) 



(27r) n q 2 — m\ — g 2 ir(q 2 ,m 2 ) 



(A186) 



This is how the function 7Tt(<?) (2. lid) appears in the present way of calculating the pressure. 



The last relation we need is the one between the coupling constants g and g = Ay A (compare the 
paragraph before (3.11a).) To order A we have 

g 2 =g 2 {l+g 2 [n'^ 2 )-7r'(m 2 )]}-\ (A.19) 

where 7r(q 2 ) is defined as n(q 2 ) in (2.12a), but with mass Hi instead of m\. From (A.19) we find for 
the derivatives of g with respect to m\ and m\ keeping g fixed: 



|| = G^KM)], = 1.2). 



(A20) 
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From this we derive easily 



-^-2 [g 2 ^T(k,ml)} = g 2 J^2 [^T(k) - 7r(ml)] 



dm\ 1 ' J dm\ 

- [k 2 -m\- g 2 it T (k, ml)] g 2 -J^tt' (m 2 ,) , 



(A.21) 



^2 [g 2 n T (k, m 2 )] = -g 2 [k 2 -m 2 - g 2 7t T (k, m 2 )] 7r"{m 2 ) (A.22) 



dm 2 



Collecting everything together we get now in a straightforward way: 



(A.23) 



where (°)p(T) mt is the interaction pressure as in (2.17) and the last term on the right-hand side is 
dm 2 , times the derivative with respect to m| of the pressure of the free photons of renormalised mass 
m2. In deriving (A. 23) we have made use of the relations 



Im Tl" (ml) = 0, 
d 



n 2 , „ (A24) 



which are valid for stable photons and electrons, that is for 

< 7772 < 4r ™i (A25) 
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